Spherically symmetric gravity coupled to a scalar field with a local Hamiltonian: the 
complete initial-boundary value problem using metric variables 
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We discuss a gauge fixing of gravity coupled to a scalar field in spherical symmetry such that the 
Hamiltonian is an integral over space of a local density. In a previous paper we had presented it 
using Ashtekar's new variables. Here we study it in metric variables. We specify completely the 
initial-boundary value problem for ingoing Gaussian pulses. 



I. INTRODUCTION 



Spherically symmetric gravity coupled to a scalar field is a rich model, where one can test scenarios of black hole 
formation, the critical phenomena discovered by Choptuik and Hawking evaporation at the quantum level. For many 
years the full quantization of the model resisted analysis, in part due to the complicated nature of the Hamiltonian 
structure of the system. Initial attempts to study the problem were done by Berger, Chitre, Nutku and Moncrief 
[l[ and further developed by Unruh 2]. The resulting complicated nature of the gauge fixed Hamiltonian led led 
Unruh to say "I present it here in the hope that someone else may be able to do something with it." More recently, 
■ Husain and Winkler and Daghigh, Kunstatter and Gegenberg Q, using Painleve-Gullstrand coordinates simplified 
somewhat Unruh's treatment. None of these efforts provided a Hamiltonian that was the spatial integral of a local 
density, leading to non-local equations of motion with the ensuing difficulty at the time of quantization. 

We recently noted that using Ashtekar's new variables the construction of a local Hamiltonian was possible. It was 
later suggested by Unruh Q and Gegenberg and Kunstatter Q that a similar construction was possible in metric 
variables. In hindsight, this is not too surprising. The key element used in our construction was that in Ashtekar's 
variables the gravitational part of the Hamiltonian constraint becomes the total derivative of a quantity with respect 
to the radial variable. It turns out that some years ago, Kuchaf [6] introduced canonical coordinates for spherically 
symmetric vacuum gravity in which one of the coordinates is the mass as function of the radius. The gravitational 
part of the Hamiltonian constraint in that case is given by the total derivative of the mass with respect to the radial 
variable. Therefore a construction similar to the one we had carried out with Ashtekar's variables can be carried out 
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■ with Kuchaf 's variables. We will detail the construction here. As most gauge fixings, only certain families of initial 

\Q ' data can be accommodated with a given choice of gauge. We set up a suitable initial-boundary value problem in the 

. gauge fixed theory and for the physically important case of Gaussian pulses. 

The organization of this paper is as follows. In section II we discuss the gauge fixing in terms of Kuchaf 's variables. 

£NJ , In section III we set up the Hamiltonian. In section IV we study the initial-boundary value problem, in particular for 

t— I ■ Gaussian pulses. We end with a discussion of possibilities for quantization. 
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II. GAUGE FIXING IN THE KUCHAR VARIABLES 

The starting point is the three-metric in spherical coordinates, 

ds 2 = A{r) 2 dr 2 + R(r) 2 (d6 2 + sin 2 Odip 2 ) , (I) 

with A and R arbitrary functions of the radial variable (and time), and their corresponding canonical momenta Pa 
and Pr. The canonical formulation in terms of these variables has been discussed by Kuchaf [f|, so we will not repeat 
it here, we refer the reader to his paper for details. The total Hamiltonian density is obtained from the Hamiltonian 
and diffeomorphism constraints, 

H T = NH + N r C r , (2) 
1_ (_A _ PaPr PjA (Rf R'RA' WR) P| R 2 

G \ 2 R 2R 2 2A A 2 A 2R 2 A 2A ' 1 ' 



C r = ^{-P' A A + P R R') + P^', (4) 



with N the lapse and N r the shift. 
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We now proceed to redefine the lapse and shift 

N R' 

N oid = ^fL, (5) 
N r oM = K ew + ^^, (6) 

and from now on we drop the "new" subscripts. The total Hamiltonian density can then be written with the 
gravitational part explicitly as a derivative with respect to the radial coordinate, 



H T = N 



1 (R'fR R P%\ R' [n2 , d4 ,^21 , PaP^' 



G \ 2A 2 2 2R 2A 2 R 2 



AR 



+ N r 



l(-P^A + P R R') + P^' 



(7) 



Let us now proceed to gauge fix. We start by setting R = r. Preserving this condition in time implies that the 
shift vanishes (this is the rescaled shift, the original shift does not vanish). One solves the diffcomorphism constraint 
to obtain P R . 

To completely fix the gauge, we need to fix another variable. With that objective in mind, it is convenient to rewrite 
the Hamiltonian as, 

with 

with Rs at the moment just a constant, later it will be identified with the Schwarzschild radius. 

The strategy for finding a gauge fixing that leads to a local Hamiltonian will be to fix the value of the quantity X. 
The resulting constraint therefore depends on the gravitational variables undifferentiated. When one preserves that 
gauge fixing in time, the lapse will be fixed by an algebraic equation rather than a differential one. This is the key 
point. If one were left as usual with a differential equation, the lapse would be an integral of the canonical variables. 
Since the Hamiltonian is an integral that involves the lapse, it becomes an integral of an integral and in that sense 
is non-local. So we choose X = f(r,t). Preservation in time of this condition determines the lapse as an algebraic 
function of <fi, P^, A, Pa. 

We proceed to solve the variable A through the gauge fixing, 

A = ~ (10) 

with 

Y=^r 2 +P 2 -rR s -2f(r,t)r (11) 
and substitute it in the total Hamiltonian, which leads to, 

_ g (g + (^) 2 r4 ) (r 2 + P l- rRs 2f(r, t)r) gP^y 

f - ^3 -2 y 12 > 

which we should solve to get Pa as a function of cj> and P^. We will see later how to do this in a compact way. 

This completes the gauge fixing. The free variables are <fr, P^. We now go to the evolution equations for those 
variables, derived before the gauge fixing, and substitute the latter. The resulting equations can be shown to be 
equivalent to those that stem from a true Hamiltonian, 

ffTruc = /(M) Gy ( PA(lp,) + y)' (13) 

with 

Y = ^r 2 +Pl{^P^)-rR s -2f{r,t)r (14) 



and in these expressions Pa should be substituted by the expressions we derived before during the gauge fixing. 
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III. OBTAINING THE TRUE HAMILTONIAN DIRECTLY 

A constructive procedure to directly obtain the above true Hamiltonian is to perform a canonical transformation 
from the variables A, Pa to a new set of variables A, P x . This should be done before the gauge fixing X = f(r, t), so 
at the moment X is function of A and -Pa given by This motivates us to consider a generating function of type 
3, F 3 (Pa, X) for which one would have that, 

a -25^Q, (15) 

and solving for A in the definition of A © this can be integrated to give 

F 3 (A,P A ) = -^l n (P A + .Z), (16) 
with Z being the non-gauge fixed version of Y, 

Z = 

We therefore have for the conjugate variable, 



Pl(^P <f ,)-rR s -2Xr. (17) 



and for A, 



gF 3 (P A ,A) = 

x dX GZ (P A + Z)' K ' 



A= dF 3 (P A ,X) = r 

dP A v/r 2 + P 2 - rR s - 2Xr 



The total Hamiltonian in terms of the new variables is 



^=N[P x + Gz{pA + z) ). (20) 



-\ 


4>, J drP to tj 
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4>, J drPtotj 
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(j), J drPtraej 



with Pa obtained by solving (jHJ with A substituted by (|T9| . 

We now proceed to gauge fix X = f(r,t). Preservation in time of this condition leads to N = f(r,t). Noting that 
cj), P,p have vanishing Poisson brackets with P x , if we write the evolution equations and substitute the gauge fixing in 
them, we have that, 

(21) 

(22) 
with 

gt - = /(M) Gr(p A (0,p,) + y) - (23) 

IV. SETTING INITIAL AND BOUNDARY DATA 

As in any gauge fixing in a complicated theory like general relativity, one does not expect one will cover all of phase 
space. The limitation here is given by the equation for Pa, which written explicitly reads, 

aPi + bPl + c = 0, (24) 

with 

» = ^( P |-r 4 (<) 2 , (25) 

b = GW (P| + r 4 (0f) - G 2 P2 (4>'f (r 2 - rR s - 2/(r, t)r) , (26) 

c = r 4 W 2 , (27) 



G(p 2 + r 4 (0')) (r 2 -rP s -2/(r,0) 



W = -f\r,t) + —± ^ . (28) 
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This will not generically yield a real value for Pa given arbitrary initial data for <fi, P^. This not only limits the initial 
data but also the boundary conditions one can give at the outer and inner boundary. So from now on we are limited 
to consider more specific situations. One has certain freedom to modify things by playing with the function f(t,r) 
that determines the gauge fixing. For different choices of f(t, r) different families of initial and boundary data will be 
acceptable as producing real values for the variables. 

A case of great interest is the study of the propagation of wave packets of scalar field on a black hole space-time. 
We will therefore concentrate ourselves on that situation. This will require specifying at spatial infinity boundary 
conditions such that the geometry is asymptotically that of Schwarzschild with no ingoing matter fields, and the inner 
boundary corresponds to a dynamical horizon with matter fields purely ingoing into it. 

We choose as initial data for the scalar field, 

= 0, r) = O ^ i 1 [r * S) , (29) 

where we are considering a Gaussian pulse and we added a factor (r — R s ) 2 /r 2 such that the field vanishes at the 
horizon initially. This makes the horizon for the initial data an isolated one. For P^ we choose what is needed to have 
a purely ingoing pulse, 



exp (~ (r r f g2 ) {r _ Rs) 2 ^ r 2 +2 _ (J 2 



P^t = 0,r)=cf> * '- v r2 °' ± - 2 ^. (30) 

We will now proceed to fix the gauge in such a way that the bi-quadratic equation (pMf has at least a pair of real 
roots. Notice that in (f2"4"| both a and c are always positive. For the quadratic equation for Pj[ have a positive root 
one needs to make the linear term negative. One possible strategy is to consider (f3"TT) and integrate it using the initial 
data we are considering, 

GlP 



(p| + r 4 (0')) (r 2 -rR s ~2f(r,t)) 



f<r, t) = ■ (3i) 

The right hand side is a bit complicated, but evaluating it numerically one sees it has the form of a Gaussian-like 
shape. One can therefore simply take for / a Gaussian that envelops the integral as a gauge choice, 

/' = Gv* (r - Rs f exp L (r = r f ^ j , (32) 

with k appropriately big for it to envelop the integrand. The integral for / can be evaluated in closed form, but its 
expression is lengthy. The form of the function is relatively simple, it is a modified step function as shown in the 
figure. In terms of / one can solve (j2"4")l for P\. The closed form expression is again lengthy. Asymptotically for large 




FIG. 1: The function providing the gauge fixing, shown for a = 0.05, G = 10 11 , Rs — 2, ro = 100. 
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r we have Pa = cor 3 / 2 with c a constant. 

This completes the determination of the initial data. We need to fix the gauge for all time. To do this, we 
consider the preservation of f(t,r) — X, which determines the lapse. We would like the lapse, at least asymptotically, 
to reproduce the usual manifestly asymptotically flat nature of the Schwarzschild space-time. This corresponds to 
N = 1 (recall that we are referring to N ncw = 1, which corresponds asymptotically to 7V id = 1/A). This results, 
asymptotically in X = X asymp . + A(r)t, where X asym p. is the asymptotic value of the expression of /, which, for 
instance, can be read off for large values of r in the figure above. The expression of A(r) is, 

A. ( 1 



A(r)= Cl r 4 exp(--(r-r o r^), (33) 



with ci a constant. 

With this form of the gauge fixing, the asymptotic form of the metric is, 



2GM , s 

9u = -1 + , (34) 



r 



9tr ^+Pl-2GMr ' (35) 

9rr = r*+Pl-2GMr' (36) 

with GM — Rg/2 + Y asymp ., which shows that the asymptotic mass is the same as that of the horizon plus the 
contribution of the scalar field. With a simple redefinition of t this yields the usual expression of the Schwarzschild 
metric in the Schwarzschild coordinates. 

Although we have not studied the evolution in detail, one can envision using a gauge with Q32I) modified to be an 
ingoing pulse and this should yield real expressions for all quantities as the pulse travels inward, at least far away 
from the black hole. 



V. DISCUSSION 



We have shown that one can gauge fix spherically symmetric gravity coupled to a scalar field in terms of the 
traditional metric variables with a Hamiltonian that is the integral of a local density in explicit form. We construct 
a family of gauge fixings that can accommodate ingoing Gaussian pulses and show that they include manifestly 
asymptotically flat coordinates. The construction of the gauge is such that it is clear that it will evolve correctly in 
time, at least for a limited amount of time. It should be noted that we have not analyzed properly the inner boundary 
condition beyond the initial slice. One presumably would like to have a dynamical horizon that increases its mass as 
the ingoing pulses progress towards the black hole, at least studying the problem classically. Quantum mechanically, 
it is less clear what one needs at the horizon, since Hawking radiation should be present. 

The question of quantization of the model in terms of this gauge implies having to deal with the quartic equation 
that generically will lead to complex values. It is therefore unclear that the constructed Hamiltonian could be 
promoted to a self-adjoint operator. It should be noted that there exist techniques Q to deal with these types of 
issues in quantization. We have recently illustrated this in a model system (8j. They are however, limited to certain 
regimes. Realistically, this type of approach is unlikely to yield insights about extreme regimes like the ones close 
to the singularity. But it may be useful in other situations, like in those in which a large black hole emits Hawking 
radiation to study, for instance, the back reaction of the weak radiation on the large black hole. 
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